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Abstract. We propose that cosmological magnetic fields generated in regions of finite spatial dimensions may manifest them- 
selves in the global dynamics of the Universe as 'dark energy'. We test our model in the context of spatially flat cosmological 
models by assuming that the Universe contains non-relativistic matter p,„ cx a^^, dark energy pq cx 0-"^"+"', and an extra fluid 
with pb 'X a"^^ that corresponds to the magnetic field. We place constraints on the main cosmological parameters of our model 
by combining the recent supernovae type la data and the differential ages of passively evolving galaxies. In particular, we find 
that the model which best reproduces the observational data when SI,,, = 0.26 is one with fig ^ 0.03, n ^ 7.68, Q.q - 0.71 and 
w ^ -0.8. 
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1. Introduction 

Recent advances in observational cosmology based on the anal- 
ysis of a multitude of high quality observational data (type 
la supernovae, hereafter SNIa, cosmic microwave background, 
large scale structure, age of globular clusters, high redshift 
galaxies), strongly indicate that we are living in a flat (Q = 1) 
accelerating Universe containing a small baryonic component, 
non-baryonic cold dark matter needed to explain the clustering 
of extragalactic sources, and an extra component with negative 
pressure, usually called 'dark energy', needed to explain the 
present accelerated expansion of the Universe (eg. Riess et al. 
1998; Perlmutter et al. 1999; Efstathiou et al. 2002; Caldwefl 
2002; Percival et al. 2002; Spergel et al. 2003; Tonry et al. 
2003; Riess et al. 2004; Tegmark et al. 2004; Corasaniti et al. 
2004). 

From a theoretical point of view, various candidates for the 
exotic dark energy have been proposed, most of them charac- 
terized by an equation of state pq - wpQ with w < -1/3 (see 
Caldwell 2002; Peebles & Ratra 2003; Corasaniti et al. 2004 
and references therein). A particular case of dark energy is the 
traditional A-model which corresponds to w = - 1 . Note that a 
redshift dependence of w is also possible but present measure- 
ments are not precise enough to allow meaningful constraints 
(eg. Dicus & Repko 2004; Wang & Mukherjee 2006). From the 
observational point of view and for a flat geometry, a variety of 
studies, and especially the SNIa data, indicate that w < -1 
(eg. Riess et al. 2004; Basilakos & Plionis 2005; Sanchez et al. 
2006; Spergel et al. 2006; Wood-Vasey et al. 2007 and refer- 
ences therein). For such a fluid the condition w < -1 is prob- 
lematic since it presents instabilities and causality problems (de 



la Macorra 2007). This may be considered as an indication that 
the dark energy interacts with another fluid, for example mag- 
netic fields. Tsagas (2001) gave an interesting perspective to 
the problem by claiming that the effect of a primordial mag- 
netic field may resemble that of dark energy through the cou- 
pling between the magnetic field and space time. 

In the present paper, we would like to investigate the po- 
tential of present day large scale magnetic fields to account for 
the effect of dark energy in the observed acceleration of the ex- 
pansion of the Universe. As we show in § 2, if the magnetic 
field is highly tangled, it cannot account for the cosmic tension 
implied by the presence of dark energy. On the other hand, we 
argue that, if the cosmic magnetic field is generated in sources 
whose overall dimensions remain unchanged during the expan- 
sion of the Universe, the stretching of this field by the expan- 
sion generates a tension that behaves as dark energy. In order 
to test our model, we introduce in § 3 an extra energy density 
term in the cosmological equations which we associate with 
the magnetic field. In § 4 we place constraints on the main pa- 
rameters of our model by performing a join likelihood analysis 
utilizing the 'gold sample' of SNIa data (Riess et al. 2007) as 
well as the ages of the passively evolving galaxies (Simon et al. 
2005). We conclude with a discussion on the possible values of 
cosmological magnetic fields in § 5. 

2. Tension in a cosmological magnetic field 

It is well known that a positive pressure in the expanding cos- 
mic fluid contributes to the deceleration of the expansion. It 
is easy to understand this when we realize that every expand- 
ing part of the Universe is pushed against the expansion by 
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its neighboring parts. As a result, each expanding part per- 
forms work against its surroundings, and thus loses kinetic 
energy. The exact opposite takes place when the expanding 
fluid feels a tension force. In that case, each expanding part is 
pulled outwards by its neighbors, and the work done on it by its 
neighbors contributes to the acceleration of its expansion (e.g. 
Harwit 1982). Dark energy is, therefore, the cosmic tension that 
accounts for the observed present accelerated expansion of the 
Universe. 

When one talks about tension, one immediately comes 
to think about magnetic fields. Magnetic field lines may be 
considered as strained ropes, with a highly anisotropic pat- 
tern of tension and pressure. We do observe magnetic fields 
up to several tens of fiG in galaxy clusters (see Carilli & 
Taylor 2002 for a review), but their origin remains a mys- 
tery. There is a general belief that cosmic magnetic fields 
are produced through some type of dynamo process that am- 
plified a weak protogalactic seed magnetic field of the or- 
der of 10^2° G (e.g. Ruzm aikin, Shukurov & Sokoloff 1988 ] 
IKulsrud ef g/. 1997l l. In this picture, the magnetic field perme- 
ates the cosmic fluid which is assumed highly conductive. In 
other words, the sources of the cosmic magnetic field are elec- 
tric currents distributed more or less isotropically throughout 
the expanding cosmic fluid. Furthermore, because of flux con- 
servation, the magnetic field scales with the expansion of the 
Universe as a^^, where a is the Universe scale factor, and there- 
fore, the magnetic energy contained in any expanding volume 
of the Universe scales as a"'. Such a magnetic field gener- 
ates a positive isotropic (on average) pressure ps - psl'i, and 
pg = l%n oc Q'"'*. The same result may be obtained after av- 
eraging out the magnetic pressure and tension terms. In other 
words, the equation of state of a highly tangled magnetic field 
is the same as that obtained for a fluid of highly relativistic 
particles, and therefore, it cannot account for the cosmic ten- 
sion implied by the presence of dark energy. The above led the 
community to conclude that isotropic tension, or equivalently 
negative pressure, is peculiar to a scalar field JZeldovich 19 86 ). 

Here, we would like to investigate the potential of a differ- 
ent scenario in which the magnetic tension does manifest itself 
in the expansion of the Universe as dark energy. We thus pro- 
pose that the sources of cosmological magnetic fields are of fi- 
nite dimensions, and that these dimensions remain unchanged 
during the expansion of the Universe. It is interesting to note 
that, galaxy clusters are the largest gravitationally bound struc- 
tures in the Universe and as such, they have decoupled from 
the background expansion. What we have in mind here is some 
mechanism that results in the generation of magnetic fields 
of order Bq around 'sources' of spatial dimensions ro (such 
as galaxies, or clusters of galaxies). Such a scenario may not 
be unreasonable. In fact, we have already proposed a physical 
picture where magnetic fields are generated without the need 
for a dynamo mechanism, through the Poynting-Robertson ef- 
fect on electrons in a highly conducting plasma around bright 
gravitating sources such as active galactic nuclei, black holes, 
neutron stars, and protostars ([Co ntopo ulos & Kazanas 1998 
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If we assume the existence of such sources of cosmic mag- 
netic fields, it is natural to further assume that the magnetic 



field around each source has a dipolar structure. In that case, 
the magnetic field drops with distance r as 



B(r)«Bo - 



(1) 



in the region influenced by a source of size ro at its center. 
We argue that the Universe may be filled with several such ex- 
panding regions of comoving size Ra{t), which contain cosmic 
magnetic sources of size ro <«c Ra. The magnetic field energy 
in each such region is equal to 



r-Ra 

Jr=ro 



BHr) 



d2 3 



(2) 



Obviously, the expansion of the Universe will gradually stretch 
each dipolar configuration described by eq. ([T]) into a monopo- 
lar configuration described by 



B(r)^Bo - 

JO 



(3) 



in the region around each sourc^H The magnetic field energy 
in the expanding region will approach the value 



B\r) 2 

47rr dr 

87r 



d2 3 

°o'o 



(4) 



which is 3 times larger than the expression in eq. (|2]l! 

We showed here that the magnetic energy Eg contained in- 
side a region of size Ra{t) increases with increasing a. This 
may be parametrized around the present epoch with a simple 
power law 



Eb = fBlrla" 



(5) 



where, / ~ 1/3, and n is a positive parameter to be determined 
from observations (see § 4). As a result, 



PB 



and 



Pb = --PB 



(6) 



(7) 



Note that in our picture, the magnetic field Bo is not pri- 
mordial, because if that were the case, by the time the Universe 
doubles its size, the magnetic field in each region of magnetic 
influence will have eff'ectively completely transformed itself 
from purely dipolar to purely (split) monopolar When that 
happens, the magnetic field contribution to the cosmological 
tension dies out. In our scenario, we expect that the physical 
mechanism responsible for the generation of the cosmic mag- 
netic field (e.g. Contopoulos & Kazanas 1998) requires a cer- 
tain number of years t in order to build a value of the order of 
Bo. 



' Actually, this is a split-monopole configuration where space is 
separated by an equatorial current sheet discontinuity across which 
the radial magnetic field changes direction. Such is the case in stellar 
magnetospheres. 
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3. Cosmological evolution 



We test our hypothesis by introducing an extra term that ac- 
counts for the magnetic field in the standard cosmological 
equations. In particular, we assume that the Universe is ho- 
mogeneous, isotropic, flat, and consists of the following three 
components denoted by subscripts 'm', 'Q' and 'B' respec- 
tively: non-relativistic matter (with zero pressure), an exotic 
fluid (dark energy), and a magnetic fluid with an equation of 
state given by eq. (|7]l. The corresponding equation of state pa- 
rameters w = PqIpq and n = -3pb/pb are assumed here 
for simplicity to be constant. Therefore, the evolution of the 
fluid densities p,„, pq and pB is given by p„, = -3Hpm, pg = 
-3(w + l)HpQ and pB - {n - 3)HpB- The scale factor of the 
Universe a{t) evolves according to the Friedmann equation: 



^2 (ay SttG^ 



(8) 



Difl'erentiating the Friedman equation and using at the same 
time the above formalism we obtain 



a AnG 

- = r-[Pm + (3w + \)pq + i-n + \)pb] 

a i 



(9) 



In this framework, we define the density parameters n„,(a), 
□^(a) and Q.b{o) as 



Q,„(a) = 



Pn 



Q.„a 



Pm +PQ+PB E^(a) 



PQ 



pm + PQ+PB 



Qea^3(i+vv) 
EHa) 



^Bia) = 



Pb 



-3+n 



Pm +PQ+PB 



EHa) 



with 



Eia) = (Q„ 



• Q.Qa 



1/2 



(10) 



Here, the Hubble parameter is given by H(a) - HoE(a), where 
Hq is the Hubble constant at the present time. Also, f2,„ + Qg + 
^B = 1. 

Note that in the context of our model, A-models correspond 
to {w,pb) = (-1,0) or (n,pQ) - (3,0), while if n = or w = 0, 
the extra fluid behaves like pressureless matter. It is interest- 
ing to mention that the interplay between the values of w and 
n could yield flat cosmological models for which there is not 
a one-to-one correspondence between the global geometry and 
the expansion of the Universe. Indeed, in a flat low-f2,„ model 
with (w,n) - (-1/3, 1) we have the same dynamics as in an 
open Universe, despite the fact that these models have a spa- 
tially flat geometry! 

In this cosmological scenario, there is an epoch which cor- 
responds to a value of c = a;, where d'/ = 0. This is called 
the inflection point. After that epoch we reach an acceleration 
phase with a > 0. Eq. (|9]l thus implies that at the inflection 
point. 



Pm,i + i3w + 1)pqj + i-n + 1)pbj = , 



(11) 



Q.m + Q.q{3w + \)af''' + Q.B{-n + \)a", = 



(12) 



Therefore, in order for the latter equation to contain roots in the 
interval of a e [0, 1], we obtain a theoretical boundary for the 
possible (n, wq) values, namely 



3wQ.Q - nD.B < -1 with w < and « > . 



(13) 



It is obvious that when ^ 0, the above constraint tends to 
the Quintessence case w < -1/3, as it should. 

4. Likelihood Analysis 

In order to constrain the cosmological parameters in our model, 
we use the 'gold' sample of 182 supernovae of Riess et 
al. 2007. In particular, we define the likelihood estimatoJl as: 
X™i^(c) (X exp[-A'2Nj^(c)/2] with: 



182 



(14) 



where M is the distance modulus M - SlogDL + 25, Z)l(z) is 
the luminosity distance Z)l(z) = (1 + z)x{z), Zi is the observed 
redshift, cr, is the observed uncertainty, and c is a vector con- 
taining the cosmological parameters that we want to fit (Riess 
et al. 2007). Note, that x is the coordinate distance related to 
the redshift through 



x(z) 



— r 

^0 Jo 



dy 
E(y) 



(15) 



We remind the reader that we work here within the frame- 
work of flat cosmology (Q^ + + ^b = 1) with non-zero 
large scale magnetic fields > 0. Furthermore, we use the 
results of the HST key project (Freedman et al. 2001) and fix 
the Hubble parameter to its nominal value //q - 72 km/s/Mpc. 
The matter density Q„, remains the most weakly constrained 
cosmological parameter. In principle, Q,„ is constrained by the 
maximum likelihood fit to the WMAP and SNIa data, but in 
the spirit of this work, we want to use measures which are com- 
pletely independent of the dark energy component. An estimate 
of without conventional priors is not an easy task in ob- 
servational cosmology. However, many authors using mainly 
large scale structure studies, have tried to put constraints to 
the Q„, parameter In particular, from the analysis of the power 
spectrum, Sanchez et al. (2006 and references therein) obtain a 
value Q,„ ^ 0.24. Moreover, Feldman et al. 2003 and Mohayee 
& TuUy 2005 analyze the peculiar velocity field in the local 
Universe and obtain the values Q,„ ^ 0.3 and ^ 0.22 respec- 
tively. In addition, Andernach et al. 2005, based on the cluster 
mass-to-light ratio, claim that £1^ Ues in the interval 0.15-0.26 
(see Schindler 2002 for a review). In the present paper, we de- 
cided to fix Q.,,, to the value 0.26. 

In this case, the vector of unknown cosmological parame- 
ters is c = (Qq, w, n). We, therefore, sample the various param- 
eters as follows: the dark energy density Q.q e [0.01,0.74] in 
steps of 0.01, the dark energy parameter w e [-4, -0.1] in steps 



Likelihoods are normalized to their maximum values. 
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2. Left panel: We show the effective equation of state pa- 
lameter as a function of the scale factor of the Universe. Right 
-panel: We present the evolution of the density parameters £1^ 
(solid), Oq (dot-dashed) and (dashed). The vertical lines 
corresponds to the present epoch. 



Fig. 1. Likelihood contours in the (Qq, w) plane for 0,„ - 0.26. 
The contours are plotted where -21nX/Xmax is equal to 2.30, 
6.16 and 11.83, corresponding to Icr, 2cr and 3cr confidence 
level. Note that the continuous and the dashed lines correspond 
to the SNIa and H{z) results respectively, while the shadowed 
area is ruled out by stellar ages. 

of 0.02, and the magnetic field scaling parameter n e [0.1, 10] 
in steps of 0.02. Doing so, the likelihood function peaks at 
- 0.5 (Qb - 0.24) with w ^ -1.32 and n ^ 0.44 (or 
wb - -0.15) which corresponds to an age of the Universe of 
12.3 Gyr The latter appears to be ruled out by stellar ages. 
Indeed, in order to put further constraints on our solutions we 
use additionally the so called age limit, given by the age of the 
oldest globular clusters in our Galaxy (^ 12.5-13 Gyr; Caputo, 
Castellani & Quatra 1988; Cayrel et al. 2001; Hansen 2002 and 
2004; Krauss 2003 and references therein). Taking into account 
the above age limit, the resulting best fit solution is: - 0.04 
(Q.Q - 0.7), w =^ -0.8 and n ^ 2.52 (or wb - -0.84), corre- 
sponding to an age of 13.1 Gyr 

In fig.l (solid lines) we present the Icr, 2cr and 3cr confi- 
dence levels in the {Q.q, w) plane by marginalizing over n. It is 
evident that Qq is degenerate with respect to w and that all the 
values in the interval < Qq < 0.74 are acceptable within the 
Icr uncertainty. Therefore, in order to put further constraints 
on Qq we additionally use measures of H(z) (see Simon et al. 
2005) from the differential ages of passively evolving galax- 
ies (hereafter H(z) data). Note that the sample contains 9 en- 
tries. Doing so, the H(z) likelihood function can be written as: 
£"(c) oc exp[-xl(c)/2] with: 



1=1 



where H(z) is the Hubble parameter (see section 3), H(z) - 
H()E(z). The dashed lines in fig. 1 represents the Icr, 2cr, and 
3cr, confidence levels in the (Q.q, w) plane. In this case, we find 
that the best fit solution is Q.q =^ 0.56 (Q^ - 0.18), w ^ -0.68, 
and « ^ 8 (with T ^ 13.2Gyr). We now join the likelihoods, 

£j™«(c) = X^'^^" X , 

and the overall function peaks at = 0.71+Q 26 (^b - 0.03), 
w = -0.80!°;i4 and n ^ 7.68!2-42 (or wb ^ -2.56). Note that, 
the corresponding age of the Universe is 13.2Gyr, while solving 
numerically eq. (12) the inflection point is at a; ^ 0.57 (or 
Zj - 0.76). Finally, due to the fact that the association of the 
extra term in the Friedman equation with the magnetic field is 
arbitrary, we may equally well consider the solution ^ 0.71 
(Qq ^ 0.03), w ^ -2.56 and n ^ 2.4 (or wb ^ -0.8). 



5. Discussion 

The above values of correspond to an average cosmic mag- 
netic field 



{B} = i^n^BPcrC") 



2.1/2 



: 650Q^^^/! A/G . 



(17) 



(16) 



where, pcr - 1.88 x lO^^^/i^ g cm"-' is the critical density of the 
Universe; h is the Hubble constant in units of 100 km/s/Mpc; 
and c is the speed of light. If ^ 0.03, (B) ^ SQfiG, whereas 
if Q.B - 0.71, (B) ^ 400yL(G. In the present work we do not dis- 
cuss what physical mechanism may account for the generation 
of such high magnetic fields. Magnetic fields of the order of 
several tens of fiG have been observed in the centers of cooling 
flow clusters. We refer the reader to Vogt & EnBlin (2003) (and 
references therein) for a detailed discussion of the techniques 
used to estimate such high values of cluster magnetic fields. 
Furthermore, a number of authors have investigated the possi- 
bility that ~ 50;uG fields may provide magnetic pressure sup- 
port in cluster atmospheres (e.g. Loeb & Mao 1994; Miralda- 
Escude & Babul 1995; Dolag & Schindler 2000; but see also 
Rudnick & Blundell 2003). We argued that dark energy (or 
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equivalently cosmic tension) manifests itself as a tension agent 
between galaxy clusters, and as such, it obviously acts in intra- 
cluster space. Moreover, there is a theoretical indication that 
strong magnetic fields lie in regions of significantly reduced 
plasma density (e.g. Gazzola et al. 2007), which would make 
observations of cosmic magnetic fields through Farady rotation 
or X-ray emission in intra-cluster cosmic voids very difficult. 
We must keep in mind that the study of cosmic magnetic fields 
on scales of clusters of galaxies is a fairly new area of research 
(see Carilli & Taylor 2002 for a review), and therefore, we can- 
not preclude future observational surprises. In particular, any 
observation of magnetic fields of the order of lOOyuG over cos- 
mological scales would give credence to our scenario. 

We would like to end this section with a short discussion on 
the equation of state of the dark energy. As we mentioned in the 
introduction, there is strong indication for an equation of state 
more complicated than the simple assumption of a constant ra- 
tio between the pressure and the energy density. We may thus 
combine our two dark energy fluids into one with an effective 
dark energy parameter 

Pq + Pb 

Weff = = . (18) 

Peff PQ+PB 

Using the evolution of pq and pB the effective dark energy pa- 
rameter as a function of time is given by 



3w - nva^"^" 
Wesia) % where v - — 

3(1 H- VQf3vv+n) Qg 



(19) 



Using our best fit parameters, Wefr ~ -0.9 for a ~ \, Weff ~ 
-0.8 in the limit a «; 1, and Wefj ~ -2.6 in the limit a » 1. 
In the left panel of fig. 2, the soUd line shows the evolution 
of the effective equation of state parameter as a function of the 
Universe scale factor. A first order Taylor expansion around the 
present epoch (Linder 2003) yields 



Weff (a) ^ -0.87 -H 0.36(1 - a) 



(20) 



In the right panel of fig. 2 we show the evolution of the 
density parameters Q,„ (solid line), Qq (dot-dashed) and 
(dashed). It is interesting that, although at the present time the 
dark energy in dominant, before the inflection point (a/ ^ 0.57) 
the Universe was matter dominated, i.e. Q.m(a <sc 1) ~ 1. 
In fact, we estimate that prior to the inflection point, D.^ < 
0.2%Q„„ which corresponds to an average cosmic magnetic 
field {B) < 3QfiG. This value was well under equipartition, and 
therefore, one may argue that, at an early enough epoch, matter 
is able to generate the cosmological magnetic fields required 
in our scenario (for example, through the Poynting-Robertson 
mechanism described in Contopoulos & Kazanas 1998). 

We conclude with a summary of the main elements of our 
scenario: 

1 . The cosmological magnetic field is generated in sources of 
characteristic size ro with characteristic value Bo- The main 
idea is that the size of these sources does not follow the 
overall expansion of the Universe. We know that the ex- 
pansion of the Universe manifests itself over length scales 
larger than the typical size of clusters of galaxies. This leads 
us to suggest that the size of our putative magnetic field 
sources is of the order of a few Mpc. 



2. Each source is associated with a region of magnetic influ- 
ence around it where the large scale field is due to the cen- 
tral source. The sources are uniformly and isotropically dis- 
tributed throughout the Universe. 

3. As the Universe expands, the magnetic field in each region 
of influence is stretched, and the total magnetic field energy 
grows. This results in the acceleration of the expansion. The 
acceleration will decrease unless new sources are continu- 
ously generated throughout the Universe. 

4. We model the effect of the magnetic field with an extra term 

oc a"^^ in the Freedman equations. The model that best 
reproduces the observational data when Q,„ - 0.26 is one 
with ^ 0.03, n ^ 7.68, Qq - 0.71 and w ^ -0.8, 
which yields an average cosmic magnetic field of the order 
of ~ 80/iG. Obviously, we may equally well consider the 
solution ^ 0.71, n 2.4, 0.03 and w ^ -2.56. 

The latter corresponds to an average cosmic magnetic field 
of the order of ~ 400yuG. 
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